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Abstract 

We use PCAC in the small region in order to calculate the Adler sum rule and the production 
of hadrons in the low energy region where resonances dominate. We find very good agreement with 
the sum rule and with the computed cross sections. We find a value C^(0) close to the Goldberger- 
Treiman prediction. The formalism is general and can be applied to other reactions shedding light 
into the dynamical transition from resonances to deep inelastic scattering. 
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I. INTRODUCTION 



Inspite of the long history of neutrino interactions with hadrons there is still interest on 
them because many cross sections, especially at low energies, are not known precisely. Anal- 
yses of experiments need the cross sections in order to interpret properties in the oscillations 
of neutrinos and separate specific quantities, like the mixing angle ^13, CP- violation, etc. 
The cross sections are also of interest on their own right because they explain properties of 
the strong interactions and the transition from the resonance to the deep inelastic region. 

In this article we consider again the production of hadrons in the small Q^-region. Matrix 
elements of the vector current have been evaluated using the conserved vector current (CVC) 
property by relating them to electroproduction data. This way the magnitude and 
dependence have been determined [1] and have been applied [2-4] successfully. For the 
vector current we use the results from [1]. The conservation of the vector current is already 
incorporated in the equations of the Appendix of [5]. 

The matrix elements of the axial current are more difficult to determine and one resorted 
to models for estimating specific values of the form factors and their Q"^ dependence. In 
addition data in the small region do not agree with some theoretical predictions [2, 6]. In 
this article we address the small region. We adapt the principle that amphtudes which 
are free of strong interaction singularities in a specific variable will generally vary smoothly 
with that variable [7]. We apply this smoothness assumption to the delta resonance by using 
the partially conserved axial vector current (PCAC) for < 5m^ ft! 0.1 GeV^ and extend 
the results to values twice as big, i. e. Q"^ up to 0.2 GeV^. This kinematic region is especially 
sensitive in charged current reactions because the three quantities Q^, and are of 
the same magnitude. As a check we calculate with the same assumptions the Adler sum 
rule and establish that it is very well satisfied for < 0.2 GeV^. For the differential cross 
section da/dQ^ we use exact kinematics by keeping the three small quantities mentioned 
above. 

We outline now the main framework and the equations that we shall use in our calcu- 
lations. The matrix element of Aj^ between a proton and the delta resonance is defined 
in Eq. (2.7) of reference [5] to which we will refer as (LP). It involves the form factors 

C^iQ^), C^{Q'^), C^{Q'^) and C^{Q'^). The last form factor is the induced pseudoscalar, 
which can be written explicitly in terms of a pion pole. When we keep the mass of the muon. 
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the leptonic current can be expanded in four polarization vectors [8, 9] . Estimates of contri- 
butions from transverse polarizations are small in the small region and will be neglected. 
The reason for this is the small Ap — )■ Trp amplitude which is primarily longitudinal [10]. 
Thus we shall concentrate on contributions from the polarizations 

= ^ (1) 

and 

eM(A = 0) = -^(|g1, 0, 0, go). (2) 

For the small region we adopt these approximations. In this case the overall tensors 
multiplying the form factors and are such that their inner product with vanishes 
(see Eq. (2.7) in reference [5]). 

We keep the two form factors and and for the pseudoscalar form factor wc write 
the pion pole explicitly. For the TrpA vertex we introduce (A | j^r | V) with j^r the pion source: 

<A++ I ^+ I p> = x/3 i^^{^)iglCt{Q^)u{v) + VS^^ (A 1 3. I p) • (3) 

q 

Taking the divergence on both sides and using PCAC we arrive at 

iq^ (A++ \A;\p) ^,{p')q'C^{Q'Mp) + (A | \ p) (4) 

q 

= V3^^(A|j.|p). (5) 
q^ — 

This relation holds for values of where PCAC is valid. It follows now that 

C^{Q')^x{p')qMp) = -U (A I J. I p) . (6) 

This is our PCAC relation which we assume to hold in the extended region of because 
there are no singularities in Q^. Obviously, we eliminated the pion pole by subtracting Eq. 
(5) from (4). 

A similar relation holds when we replace the delta resonance by a general final state 
In this cases 

{x-^ I a; I p> = + v^-^ <^++ I A I P) (7) 

q in^ 

whose divergence gives the relation 

q^n^^-tU{X++\j^\p). (8) 



Here TZfj_ is the remaining amplitude beyond the pion pole. 

The reader should note that this relation holds for any final state mx 7^ fUp and also 
when we replace the proton in the target by a neutron. We need this relation for estimating 
the reaction A{W^n — )■ X^) = A{W~p — )■ X^) in the sum rule and also for contributions 
uptoW < 1.6 GeV. 

Another remarkable property of neutrino interactions is that for small values of relative 
to the neutrino energy Ei, and energy transfer 1/ the leptonic current can be approximated 
with the two polarizations given in Eq. (1) and (2). The scalar polarization gives a small 
term proportional to which we keep. The polarization with zero helicity annihilates the 
pion pole and gives the dominant contribution for the remainder as follows: 

e'^iX = 0) (A++ \A,\p) = e^{\ = Q)iVM,{p')Ct{Q')u{p) (9) 

« iV?,-^^,{p')Ci{Q')u{p) + O (10) 

A{t:+p A++) (11) 

with = 0.093 GeV. We emphasize that we first eliminate the pion pole and then approxi- 
mate the helicity zero polarization by q^/ \fQ^. This method leaves out terms of O [Q"^ /v'^) 
and justifies the apphcation of PCAC. Thus the polarization with zero hehcity takes the 
divergence of the axial current. Using this result we calculate in the following sections the 
differential cross section and the Adler sum rule as a function of Q^. 

In previous articles the axial contribution was accounted for by introducing form factors 
and estimating the dependence from neutrino data. The Goldberger-Trciman relation 
that was obtained earlier [1, 11] is also inherent in Eq. (6) when we substitute bare couplings 
for the matrix elements. What is different now is the replacement of the left-hand side of 
Eq. (6) by 7r+p — >■ or 7r+p — >■ data. As it will become clear later on, PCAC with 
the smoothness assumption provides an estimate of the form factor for < 0.2 GeV^. The 
handling of the vector-axial interference (the structure function Wz{Q'^, u)) will be explained 
in section 3. 

In an extension of the model some authors [3, 4] , in order to account for the non-resonant 
background, use additional diagrams generated by an effective Lagrangian. In this article we 
compute the contribution of the axial current using PCAC and for the vector contribution 
we use the results from (LP). The production of pions through the exchange of mesons and 
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Regge trajectories in the t-channel has been also worked out [12]. Our PCAC results in 
Eqn. (6)- (11) sum up the contributions of these exchanges. 



II. GENERAL FORMULAS 



From the results of the previous section we determine the contribution of the axial current 
alone to the cross section 

2" 



dQMi/ 



27r 47r El 



+ Lii 



mz 



(12) 

The matrix elements Lqq, Liq, Lu were introduced and elaborated in references [8, 9]. The 
elements Liq and Lu are proportional to and with an additional factor of or 
become too small and will be neglected. The cross section a{'K'^p — )■ X'^^) is the production 
of hadrons at the energy v. It is evident now that integrating over v in order to obtain 
dcr/d(5^ we trace the delta resonance peak and this influences the turning over of the cross 
section at low Q^- 

For evaluating the Adler sum rule we need the structure function W2{Q'^,v). The identi- 
fication is more evident when we consider the dominant contribution in our Eq. (12) in the 
limit Ey — V ^ V, The leptonic density matrix element squared was computed in 

[8] as 



^00 



[Q\2E, -u)- uml 



212 



Q2(g2 + ^2) 

which in the high energy limit becomes 

2g2 



(13) 



-^00 — >■ 



m, 



4E^E'-{Q' + ml)--^u' 



g2 + i/2 

By comparing our Eq. (12) with the W2{Q^, i') in Eq. (2.10) of (LP) we find the relation 



(14) 



TT Q2 + ^^2 

with the superscipt A denoting the axial contribution alone 



(15) 
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III. TESTING THE ABLER SUM RULE 



A basic relation in particle physics is the Adler sum rule for neutrinos [13]. It follows 
from current algebra and holds for each value of Q^. For the axial current it takes the form 

POO 

[gA{Q')T + / du [Wl^iQ', v) - Wl^{Q\ v)\ = 1 (16) 

with gAiQ"^) the form factor for the vertex | ^+ | n). 

With the functional form for W^IQ"^, u) given in Eq. (15) we obtain 

[9a{Q')] ' + ^ r d^T^^ k"n^) - = 1- (17) 

The functional form in Eq. (17) follows from the general formalism. It is a very convenient 
form because at = it reduces to the Adler- Weisberger relation [14, 15]. 

We tested the relation by using experimental data for the pion-nucleon scattering and 

sAQ') = (18) 

with Ma = 1.0 GeV. 

The data we use are from the particle data group [16] and from Barashenkov et al. [17]. 
Both sources give compilations from a large number of experiments. The data gives both the 
elastic and total cross sections for the reations Tr^p and 7r~p which is equal to 7r+n. The values 
of the cross section did not change in the last few decades. In the region of integration the 
inelastic channel for TC^p is negligible. For TT~p it is substantial which indicates the presence 
of non-resonant background. For the sum rule we use the total cross sections. In the next 
section we shall use data for the Tr+p channel for extracting C^, where as mentioned the 
difference between elastic and total cross section is neghgible. 

In order to give some typical value we tabulated in Tab. I and II values of the cross 
sections for two different values of = 0.010 GeV^ and = 0.150 GeV^. The cross 
sections correspond to the values of W. In the last column we give values for the variable 
^ = '^E^Q2 which enters Eq. (12) as an overall factor. The cross section is the product of A 
with the pion-nucleon cross section which indicates that the overwhelming contribution us 
close to the delta resonance peak. 

The results of the numerical integration are shown in Fig. 1 with the resonance contri- 
bution corresponding to the integration of the total cross sections. The integral has been 
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V [GcV] 


W [GcV] 


(j(7r+p) [mb] 


aiix p) [mb] 


[GeV] 


0.20 


1.118 


16 


12 


0.329 


0.25 


1.159 


77 


30 


0.240 


0.30 


1.199 


189 


67 


0.171 


0.35 


1.238 


175 


63 


0.119 


0.40 


1.275 


95 


37 


0.079 


0.45 


1.311 


60 


28 


0.047 


0.50 


1.347 


42 


26 


0.022 


0.55 


1.381 


31 


28 


0.002 



TABLE I: Values for = 1 GeV and = q.OIO GeV. 



truncated aX v = 1.6 GeV. The sum of QE and RES should saturate the sum rule. The dif- 
ference 1.0 - (QE+RES) corresponds to contributions from higher energies, where multipion 
production is important. This is indicated in the figure as HE. The HE contribution is very 
small for ~ 0.0 and increases monotonically indicating that the multipion contribution 
is more important for larger values of Q^. Bodek and Yang, after analyzing data in deep 
inelastic scattering [18], studied the transition to the resonance region and suggested an 
interpolation, which agrees rather well with our estimated value from the sum rule. The 
sum of the three contributions produces the value of one to 10 % or better. The errors in 
the hadronic cross sections are less than 3 %. Thus the confirmation of the sum rule is an 
indication for the validity of PC AC. 

IV. THE DIFFERENTIAL CROSS SECTION 

For the calculation of the various terms we apply different methods because we choose 
for each current the method that is most reliable. For the axial current we use the PCAC 
result in Eq. (12). For the vector current we use the form factors from reference [1]. The 
vector-axial interference term is obtained by using formula (A3) in LP with the vector form 
factors from [1, 5] and the axial form factor C^{Q'^) extracted later on in this section. We 
elaborate on each of these terms. 

1. For the axial contribution we use the data for pion-nucleon cross sections [16]. We 



7 



V [GcV] 


W [GcV] 


(j(7r+p) [mb] 


aiix p) [mb] 


[GeV] 


0.20 


1.053 


16 


12 


0.109 


0.25 


1.097 


77 


30 


0.113 


0.30 


1.139 


189 


67 


0.111 


0.35 


1.180 


175 


63 


0.104 


0.40 


1.219 


95 


37 


0.095 


0.45 


1.257 


60 


28 


0.085 


0.50 


1.294 


42 


26 


0.074 


0.55 


1.330 


31 


28 


0.063 


0.60 


1.364 


24 


28 


0.052 


0.65 


1.398 


17 


34 


0.042 



TABLE II: Values for = IGeV and = 0.150 GeV. 
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[GeV^] 

FIG. 1: Solid: Adlcr sum; dashed: form factor; dot-dashed: integral over resonances; dotted: high 
energy part (Bodek-Yang interpolation). 

concentrate on the reaction v^p — ?> fi^pn^ where the non-resonant background is 
negligible so that we select X^~^ ~ A^"*". We integrate the cross section over v tracing 
the product Loo(Q^, Ei,)aTjp{v) and we obtain the curve denoted as PCAC in Fig. 2. 
Our method shows that the turn over in is in part a reflection of the peak in the 
invariant mass of the resonance. In the same flgure we show the contributions from 
{C^Y, the interference term and the remaining form factors (rest) which are smaller. 

We also determine the axial form factor C^iQ"^) by computing the integral 
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FIG. 2: Contributions to the total cross section at E,^ = 1 GeV with no cuts. 

■fi'th^'^^ ^2^^^^^' ^^^^ numerically using Eq. (15) and as a second method we use Eq. 
(A2) of (LP) where only and are kept. The two results contain axial-current 
contributions and by comparison we extract values for C^{Q^) as shown in Fig. 3. 
For data we used the reaction vr+p — >■ A++ where the non-resonant background is 
smallest. It is reassuring that the value C^(0) = 1.08 is close to 1.20 predicted by 
the Goldberger-Treiman relation. In the same figure we plotted C^{Q'^) from the 
parametrization in (LP), noting that the two results are close to each other. Several 
recent articles calculated C^(0) by fitting the experimental data [3, 4, 19-23] with 
their values varying from 0.87 up to 1.20. Models with a resonant background [3, 4] 
prefer the power value, while the other articles [19-23] prefer values closer to 1.20. 
The reasons for the differences is the treatment of the non-resonant background and 
the exact kinematics at small Q"^. 

2. For the vector contribution we use the formulas from [1] 

C3W) = |^ , C-(Q') = =^^ (19) 

QW)--^^ and D.-(l + ^X (20) 

with My = 0.84 GeV, which have been extracted from electroproduction data. The 
result is also shown in Fig. 2. 

3. For the vector-axial interference W3{Q^,u) we use the form factor C^{Q'^) extracted 
through PCAC, — —'^/aC^ and the vector form factors just described. Then we 
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1.4 
1.2 
1.0 
0.8 
0.6 
0.4 

0.00 0.05 0.10 0.15 0.20 

^ [GeV^] 

FIG. 3: Extracted values for proton scattering (solid) and the phenomenological model 
C^{Q'^) = (^i_^_Qyj^^'2 -)2 i+2q2/m^ with Ma = 1.05 GeV (dashed). Throughout our calculation we 
use the extracted values. 

use the formulas from (LP) and calculate their contributions to the cross section. The 
three contributions are shown separately in Fig. 2. 

We note that the three contributions are comparable and are important for determining the 
dependence of da/dQ"^. 

V. COMPARISONS 

For comparisons with data we shall account for two bubble chamber experiments [24, 25] 
and the recent results from MiniBooNE [26]. The results of our calculation using the Argonne 
experiment (ANL) arc shown in Fig. 4. We weighted the theoretical curve with the ANL flux 
and limited W < 1.4 GeV. A similar comparison was performed for the Brookhaven (BNL) 
experiment shown in Fig. 5. The experimental data is given in terms of event rates [25]. For 
comparison with absolute cross sections we use the transformation coefficient obtained in [4] . 
In both cases the curves are close to the experimental points. We emphasize that we did 
not include the non-resonant background, which has been estimated in electroproduction to 
be 10 % [27, 28]. 

The MiniBooNE group also reported data for a CH2 target and an absolute normaliza- 
tion [26]. For the nuclear targets we use the incoherent sum of scatterings on neutrons and 
protons. We use for each channel the same formulas and the appropriate hadronic cross 
sections. For this we calculate also the cross section on neutron targets using hadronic data. 
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■ ^ rest 

0.00 0.05 0.10 0.15 0.20 

[GeV^] 

FIG. 4: Contributions to the total cross section for the ANL flux with W < 1.4 GeV. 




0.00 0.05 0.10 0.15 0.20 

[GeV^] 

FIG. 5: Contributions to the total cross section for the BNL flux with W < 1.4 GeV. 

For the ratio of neutrino scattering we found the ratio Onl^^p shown in Fig. 6. We weighted 
the cross sections with the neutrino spectrum and show the results in Fig. 7 by taking 8 
proton and 6 neutron cross sections. Again we did not consider explicitly the non-resonant 

background or any rescattering (FSI) that may occur in the carbon target. Note that the 
results in Fig. 7 are for the entire molecular target which makes the cross section very big. 

Finally we consider a very high energy reaction from the FNAL 15ft bubble chamber [29] 
and the results are shown in Fig. 8. 

VI. SUMMARY AND OUTLOOK 

We applied the PCAC relation to the small region and calculated the contribution 
of the axial current. An interesting result is the expression for Wi{Q'^ , v) in terms of pion 
nucleon scattering cross sections, Eq. (15). We use the new expression to compute the 
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0.8 - 



V 0.6 
0.4 



0.2 - 



0.0' — 
0.00 



0.05 



0.10 
2^ [GeV^] 



0.15 



0.20 



FIG. 6: Ratio a^/a^ for the MiniBooNE flux. 




0.20 



FIG. 7: Contributions to the total cross section for the MiniBooNE flux with W < 1.4 GeV. 

Coherent sum of proton and neutron cross section for CH2. Neutron cross section are scaled 
proton cross section by ratio of direct axial contribution ratio of proton and neutron data. 




0.20 



FIG. 8: Contributions to the total cross section for the FNAL-15ft experiment with W < 1.4 GeV. 
Neutrino energy flxed at 25 GeV. 
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Adler sum rule and showed that it is saturated almost completely by data in the delta 
energy region, leaving a smaller contribution from higher energies, which becomes more 
important as increases. The higher energy contribution suggests an interpolation given 
in footnote [18]. 

We also studied extensively a general formula for the differential cross section which 
includes the mass of the charged lepton exactly. Our estimates of the cross section for charged 
pion production by neutrinos shows that contributions from vector, axial and interference 
terms are comparable. The addition of all three determines the magnitude, position and 
shape of the peak in the differential cross section at < 0.20 GeV^. Encouraged from 
the results wc computed the differential cross sections for the Argonne, Brookhaven and 
MiniBooNE experiments. The comparisons in Fig. 4-7 are very good. Besides the low 
energies a calculation for — 25 GeV produces the measured cross section. We did not 
include the non- resonant background which can be the subject for further studies [4]. For 
this reason we concentrated on the A''""'' production where the background is the smallest. 
For medium heavy nuclei rescattering corrections will also be required. 

Our method allows to extract the axial form factor C^{Q'^) whose value at = is 
consistent with the Goldberger-Treiman relation. This form factor has been the subject of 
several articles [3, 4, 19-23] with its value varying from 0.87 to 1.20. We feel that a precise 
treatment of the various contributions and of mass effects are essential. 

Our approach can be extended and combined with experimental data in order to predict 
cross sections at other energies, for reactions with antineutrinos and also for neutral currents. 
The kinematic region we considered is the place where coherent scattering on nuclei also 
occurs. Our cross section should be useful in subtracting the resonant contribution leaving as 
a remainder the coherent production. The PCAC relation and the cross section in Eq. (12) 
are general and can be applied to other resonances. 

Finally, the rapid saturation of the sum rule and estimates of the cross sections hint 
how the transition from resonances to the deep inelastic region is attained dynamically. 
The picture that emerges for da/dQ^ includes a constant contribution from the algebra of 
commutators plus a term varying with from the form factor and low mass resonances 
which fade away as increases, leaving space for multiple final states. Our results support 
and extent the results obtained earlier [30]. 
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